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Abstract 

We compute the 0(a s ) corrections to the Wilson coefficients of the dimension five 
operators emerging from the Operator Product Expansion of inclusive radiative B 
decays. We discuss the impact of the resulting 0(a s AQ GD /m%) corrections on the 
extraction of m& and ^ from the moments of the photon spectrum. 
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1 Introduction 



The inclusive radiative decays of the B meson play a central role in the search for new physics. 
While the total rate of B — > X s 7 is sensitive to new physics in flavor- changing transitions, 
its photon spectrum is almost completely determined by Standard Model physics and can 
be employed to extract information on the B meson structure and on the b quark mass, 
see for instance [U|2]. The latter are useful in a number of B physics applications, like the 
determination of \V u b\ and of |V^,|, which are both important inputs for the determination 
of the CKM unitarity triangle and the study of CP violation in the Standard Model. 

After integrating out the heavy degrees of freedom at the electroweak scale and evolving 
the resulting weak Hamiltonian to the b quark mass scale, inclusive radiative and semilep- 
tonic B decays are well described by an Operator Product Expansion (OPE) in inverse 
powers of the b quark mass. In this way one can factorize the long distance dynamics into 
the matrix elements of a few local operators [Hill]. Since the Wilson coefficients of these 
operators are perturbative, and the matrix elements of the local operators parameterize the 
non-perturbative physics, the radiative and semileptonic rates and the moments of their dis- 
tributions are double series in a s and A/m^, with A being the QCD scale. The lowest order 
of this expansion corresponds to the decay of a free b quark and linear 0(A/rrib) corrections 
are absent. The relevant parameters are therefore the heavy quark mass mj (and possibly 
m c for the charmed decays), a s , and the matrix elements of the local operators: \i\ and (Aq 
at 0(l/ml), p 3 D and p\ s at 0{l/m\). In the case of b — )■ 57 transitions, the first and second 
central moments of the photon spectrum, (E y ) and (E 2 — (E^) 2 ), are proportional to the 
b quark mass and to /i 2 , respectively, up to subleading corrections in the l/m b expansion, 
and contribute in an important way to the global fits for the extraction of |V^,|, m& and the 
other OPE parameters (see [5] for recent results). 

The precision of B — > A s 7 calculations in this framework is known to be limited in two ways: 
i) The dominant contribution to the radiative b decay is associated to the electromagnetic 
dipole operator 7 = (a e m/47r) mbSLcr^F^bn, with additional operators appearing at 0(a s ). 
It turns out that a local OPE can be written down for the O7 contribution only, while 
the other operators are expected to induce unknown 0(a s A/m&) contributions, see e.g. 
[6]; ii) Even for the dominant Oj component, the experimental cuts employed to isolate 
the signal introduce a sensitivity to the Fermi motion of the b quark inside the B meson 
and tend to disrupt the OPE. One can still resum the higher order terms into a non-local 
distribution function |7j and since the lowest integer moments of this function are known, one 
can parameterize it assuming different functional forms [8] , although alternative approaches 
are also possible [I]. 

The perturbative corrections to the leading O7 contribution are presently known to com- 
plete Next-to-Next-to-Leading Order (NNLO), i.e. 0(a 2 ), while the NNLO contributions 
involving other operators, quite important for the rate, are not yet complete [H], see JTU] 
for recent updates. As for the non-perturbative corrections to the O7 contribution, they 
are known through 0(1 /ml) [Tf]IT2"]. Among the non-perturbative corrections associated to 
operators other than O7, only one class is known: they are of 0(l/m 2 ) and 0(l/(m 2 mb)) 
and numerically small [12]. Overall, the known power corrections modify the total B — > X s 7 
rate by just ~ 3%, but they are essential in the calculation of the moments. 

In this paper we present the first calculation of the 0(a s ) corrections to the A 2 jm\ terms in 
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B — > X s 7 . The 0(a s ) perturbative corrections to the fi^ coefficient in the total rate are fixed 
by Lorentz invariance [3J. Their contribution to the photon moments has been computed in 
the context of a multiscale OPE [13] which applies however to the end-point region only. We 
have computed the relevant Wilson coefficients at 0(a s ) by expanding off-shell amputated 
Green functions around the b quark mass shell, and by matching them onto local operators 
in Heavy Quark Effective Theory (HQET). Our results allow for an improved analysis of the 
radiative moments. In particular, the inclusion of the 0(a s ) perturbative corrections to the 
variance of the spectrum, permits the extraction of fi^ at Next-to-Leading-Order (NLO). 

The outline of this paper is as follows. In section 2 we recalculate the Wilson coefficients 
of the matrix elements of dimension five operators, /ij; and /j/q, at leading order in a s and 
introduce our notation. Section 3 is devoted to the calculation of these Wilson coefficients at 
0(a s ) . It also contains our final results and a first numerical estimate of their importance. 
In section 4 we summarize and conclude. 



2 Leading Order Coefficients 

At low energy the b — > S7 transition is governed by the effective Lagrangian 

C cS = £qcdxqed(«, d, s, c, b) + ^V tb V* Cf(fi)0 7 + ..., (2.1) 

which includes, apart from the QED and QCD interactions of the light quarks u, d, s, c and 
b, the renormalization scale dependent effective Wilson coefficient C7 and the corresponding 
electromagnetic operator O7. The latter mediates the b — > 37 transition and is given by 

7 = -^m b Qi)(8o*'P R b)F ia , , (2.2) 

where m^n) is the running b quark mass in the MS-scheme, = i [7^,7^/2 and Pr^l = 
(1 ± 75)/2. The ellipses in equation (12. ip denote contributions from operators that are not 
relevant for our calculation. 

Using this effective Lagrangian we can derive that part of the differential decay rate that 
is induced by the self-interference of the electromagnetic dipole operator, 

iTll[ B _> x. 7) = I V tb V>f\Cf Ml 2 ^ «W~* . (2.3) 

Here, mg is the mass of the B meson, q the momentum of the photon, the hadronic 

tensor, 

d 3 Px 



= E / (4§k (27r)V4)(ps ~ px " q) 



x (B(p B )\XO)a^P L s(0)\X s ( Px )}(X s (p x )\m^PRb(0)\B(p B )} , (2.4) 
and pv va fi the photon tensor, 

p^aP = J2 (0|FH 7 (g,A))(7(g,A)|F^|0> 



A=±l 
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Figure 1: The imaginary parts of these tree-level diagrams contribute to the Wilson coeffi- 
cients of the operators with dimension 3, 4 and 5. 



= q a q v g^ - q a q^g vP + q^q P g va - q t Y ' g m ■ (2.5) 
The hadronic tensor W^ ya p itself can be rewritten as 

Wpvce = 2Im (i j d 4 xe- iq ^B(p B )\T{b(x)a^P L s(x)m^PRb(0)} \B(p B )) 



(2.6) 

which is useful because the time-ordered product can be expanded into a series of local 
operators that are suppressed by powers of the b quark mass. Hence we can write 

WnvapP^ " 13 = — 167Tm b fcdim3Cdim3 H Cdim40dim4 H o C dim50dim5 + •••)' (^-7) 

where Odimn is an operator of dimension n that contains n — 3 derivatives, and 



.(0) , ^£ „(1) 
-dimn ^ 



+ — c (1) + f2 



is the corresponding Wilson coefficient that can be determined in perturbation theory. In 
the reminder of this section we will review the calculation of c^ mn for n = 3, 4, 5. 

The two Feynman diagrams contributing to these Wilson coefficients are depicted in Fig. HJ 
The momentum of the incoming b quark is m^v + k, where v is the velocity of the B meson, 
and A; is a residual momentum that accounts for the interaction of the almost on-shell b 
quark with the light degrees of freedom in the B meson. The components of this residual 
momentum are of 0(A), i.e. much smaller than the b quark mass. The same is true for the 
components of the momentum r of the radiated soft gluon that is present in the Feynman 
diagram on the right-hand side. To achieve the anticipated OPE (12.71) we perform a Taylor 
expansion of the amputated Green functions corresponding to the two tree-level diagrams 
up to 0(k 2 ) for the left diagram and up to 0(k) and 0(r) for the right diagram. This is 
equivalent to an expansion in inverse powers of mj, or more precisely, to an expansion in 
(iD p — nifyV fj) j nib, where _D M = d^ + iggG^T 11 . In the algebraic manipulations we refrain from 
using any on-shell relations for the b quark, that is we do not impose any restriction on the 
residual momenta. The only on-shell condition we apply is that of the B meson, namely 
v 2 = 1. After extracting the imaginary part of the amputated Green functions we replace the 
residual momenta k and r with derivatives acting on the b and the gluon fields, interpreting 
the latter as parts of a covariant derivative. We end up with expressions for the two tree- 
level diagrams in terms of local operators of the form bF(v, q, D)b, where F(v, q, D) is a 
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matrix- valued function of v^, and in spinor space. As far as the s-quark is concerned, 
we will use m s = 0. 

At this point it is convenient to introduce HQET [H] which incorporates the expansion in 
{iDp — m^v^jm^ in a natural manner. Its Lagrangian is given by [T5|[T6] 



£hqet = ib v v ■ Db v + ^—b v (iD ± ) 2 b v - a(fi)^-b v cx^ u G^b v + O 



mt 



2 / ' 



(2.9) 



where = D» - v" v • D and [T7] 



a(fi) = 1 + 



C, -C A ( 1 + ln-^ 

m b 



a s 

7 ± + --- • 

47T 



(2.10) 



The numerical values of the color factors are Cp = 4/3 and Ca = 3. The relation between 
the b quark fields in QCD and in HQET reads 



b{x) 



ipi 



2m> 



1 + )b v (x) + 0[ -j 



1 



mi 



2 I ' 



(2.11) 



at tree-level as well as at Of a. 



Using ( 12. lip , together with fb v = b v , in all local operators of dimension 4 and 5, we end 
up with the following set of operators, 



Of = b v iD% v 



0» v = b v l -{iD\iD"}b v , 
Of = b v ^G afM a cr au T a b v 



(2.12) 



where g s G° p . v T a = — i[Dfj,, D v ]. We also find operators that include a 75. However, we will 
eventually calculate matrix elements of operators between B meson states. From parity 
considerations it follows then that only the operators given in ( I2.12p give non-vanishing 
contributions. Hence it is not necessary to calculate the Wilson coefficients of those operators 
that include a 75. We should also mention that up to this point of the calculation we did 
not use the equation of motion for the b v field that follows from the HQET Lagrangian. 

For the Wilson coefficients of the operators given in (12.121) we find, in d = 4 — 2e space-time 
dimensions, 



,(o) 

JO) 



(l-e)q fi S(l-z), 

(1 - e) [ (2u„ - q,) 5(1 -z) + (v, - q,) 5'(1 - z) 



c 2°l = (1 - e ) 2 idiw + v^q„) 5(1 - z) + (g^ v + 6v,j,v v - 6v^q u ) 5'(1 - z) 



+ 2 (v^Vy - 2v li q v + q^) 5"(1 - z) 



"3^j,y 



1 + e 
*»vflv ^—9ti, 



6(1 -z)+ (3 + e)q^q u -2v^q u + 



1 + e 



-9^ 



S'(l-z), (2.13) 



where z = 2 v ■ q = 2E 1 /m b and q = q/m b . 
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In the last step we have to calculate the forward matrix elements of the four operators given 
in ( I2.12p between B meson states. Since we expressed the leading operator O b in terms of 
the b quark fields of QCD we know its matrix element exactly, 



(B{p B )\0£\B(p B )) = 2m B vr 



(2.14) 



The evaluation of the matrix elements of the other operators involves the equation of motion 
of the effective theory, and leads to two additional matrix elements 



Ai 



2m 



B 



-(B(v)\b v {iD) 2 b v \B(v)) , A 2 



■ -^(B(v)\b v ^-G^b v \B(v)) 



(2.15) 



The velocity dependent B meson states used here are related to the momentum dependent 
ones introduced in (12.41) by \B(p B )) = ^m B \B{v)) + 0(l/m&). The power correction to this 
relation is irrelevant for our calculation. While Ai and A 2 are defined in the asymptotic HQET 
regime, in practical applications one deals with /j, 2 = — Ai+0(l/mb) and /i^ = 3A 2 +0(l/m6), 
defined in terms of the m^-finite QCD states appearing in (12. 4p . 

Now we are in a position to calculate the contribution to the differential decay rate due to 
the electromagnetic dipole operator only. We obtain 



dT 



77 



dz 
where 



,(o) 

77 



+ C 



(o) Ai 



.(0) 



Al 2m? + ° Aa 2ml 



A 2 (aO oi s (/i) f (i) (i) Ai 



47T 



+ c 



A 2 (/^) 



Al 2ml + ° X2 2ml 



(2.16) 



,(o) _ G 2 F a em m 2 b (n)m 3 b 



V tb V t *f\Cf(»)\ 



77 



32tt 4 



(2.17) 



and 



J0) 

"A 2 



6(1 



cf ) = 5{l-z)-5\l-z)--5"{l 



-9 5(1- z) -3 5'(1- z) 



(2.18) 



This is in agreement with the well-known results given in [IT] . The calculation of the Wilson 
coefficients and is the subject of the next section. 



3 Next-to-Leading Order Coefficients 

In order to determine the Wilson coefficients Cg , c^ 1 and c^ 1 we calculate the amputated 
Green functions corresponding to the Feynman diagrams shown in Fig. [2j The momentum 
assignments of the external lines are exactly the same as in the tree-level calculation per- 
formed in the last section. For the Green functions that contain a radiated soft gluon we 
apply the background-field formalism [19]. Furthermore, we work in the general i?^-gauge 
for the gluon propagator. Again we perform a Taylor expansion up to 0(k 2 ) for the diagrams 
without a radiated soft gluon and up to O(k) and 0(r) for the diagrams with a radiated soft 
gluon, and refrain from using any on-shell relations for the b quark. We apply integration- 
by-parts techniques [20] to reduce all integrals to a few so-called master integrals and solve 



6 



Figure 2: One-loop diagrams contributing to the Wilson coefficients of the operators with 
dimension 3, 4 and 5. The momenta assignments of the external lines are the same as 
in figure \T\ Sixteen additional diagrams with a gluon radiated off an internal line can be 
obtained from the diagrams shown here. 



the latter analytically. Ultraviolet as well as infrared divergences are handled by dimensional 
regularization. The ultraviolet divergences can be removed by appropriate renormalization. 
For the self-mixing of the operator 7 we use the MS-scheme, 

yMSyUS _ i i Cf «s / q . x 

and for the field renormalization constant of the b quark we apply the on-shell scheme, 

Z° s = l-C F (- + 4 + 6\n-^)^ + .... (3.2) 
\ e m bJ 4tt 

As far as the renormalization of the background gluon field G a is concerned, we only have 
to remember that g s G a ^ does not get renormalized. Finally, we extract the imaginary part 
of the amputated Green functions, replace the residual momenta k and r with covariant 
derivatives iD^ — m b v^ and the b quark spinors with b v spinors via ( 12. lip , and calculate the 
forward matrix elements of all operators between B meson states. Our result is then a linear 
combination of the three matrix elements introduced in ( 12.141) and ( 12. 15ft with coefficients 
that are ultraviolet finite but still contain infrared divergences, 

fo ( z , €, A — J v » + hi ( z i €, H, — J tt~ + h a (z, €, H, — J 7T~ ■ ( 3 - 3 ) 

We have also made explicit the dependence on the gauge parameter £. The infrared diver- 
gences are removed in the matching procedure. 
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We now turn our attention to the right-hand side of the matching equation (12. 7p . At the 
one-loop level it schematically looks like 



,™r 3 

n=3 



c dimn(°dimn)l-loop + ( ^ + 8Z, 



.(1) 



c (0) 

dim n (dim n 



(3.4) 



where (Odimn)tree and (Odi mn )i-ioo P denote the tree-level and one-loop matrix elements of the 
operator Odimn between B meson states, respectively, and Zdimn = 1 + 5-^dimn collects the 
Z-factors to render this expression ultraviolet finite. In the case at hand we only have to 
consider the one-loop matrix elements of the operators given in (12.121) since these are the only 
ones that have non-vanishing Wilson coefficients at the tree-level. The same holds for the 
tree-level matrix elements that are multiplied by Z-factors. Because of the Taylor expansion 
in the momenta k and r the one-loop matrix elements of the operators 0% and O3 \ vanish 
in dimensional regularization and we only need to compute the one-loop matrix elements of 
O^. The operators that include a 75 can again be discarded due to parity considerations. 
For the renormalization constants, we use the on-shell scheme for the b and b v spinors, and 
the MS scheme for the operator renormalization, 



b J 
fjbi bare 



[c ltl O?r e = c lM Of , [c 3 , u 3 
A simple one-loop calculation yields 



2 J 

jx ui bare 



yOS yMS^ua/3 n 
6u A, : „ C2nv u 2al3 1 



b v kin 

•OS 7 MS,/M/a ( 8 n 

b v ^chromo c 3Ati' Ly 3a/3 



(3.5) 



yMS. fiua/5 
^kin 



yMS, fiua/5 
chromo 



-C F - — - (g^ - 2v"v v ) v a v^ — + 
e Att 



— (g» a - V » V a ) g»P — : ... 



4tt 



(3.6) 



The Feynman gauge is obtained by setting £ = 1. The on-shell renormalization constant 
can be set equal to 1 since the 0(a s )-corrections to the selfenergy of the b v field depends 
on no other scale than the renormalization scale /x. 

Requiring equality between (I3.3P and (13.41) . and solving for c^ mn , we obtain infrared finite 
and gauge independent expressions. Writing the coefficients as follows, 



A 2 



(MO 



C, (c£ F) + Ar 



(1,F) 

A 2 



C A c 



,(1,A) 
'A 2 



+ Ac 



(1,A) 

A 2 



(3.7) 



our final results are given by 
4tt 2 



,(1,F) 



5 + 



3 

ln(l - z 
1 - z 



5(1-*) -7 



1 



l-z 
2 



7 + z-2z 2 -2(l + z) ln(l-z) 



4ci 0) lnA 



(3- 



,(1,F) 



--(15 + 27T 2 ) 5(1 -z)- -(3- vr 2 ) *'(!-*) 
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7 4vr : 



6 9 
1 

1-z 



2 


1 


5 


1 


3 


[(1 -z)\ 


+ ~3 


[(1-^) 2 J 



+ 



ln(l - z) 



ln(l - z) 



-4 



ln(l - z) 



,(1,F) 

-A 2 



41 
-9 
+20 

2(7 



1 - ^ 

207T 2 



+ 10 + - z - - (4 + 3z) ln(l - 2) - 4cSf In £- 



1 



3 3 
^5(1 - z) + 4 (1 + 7T 2 ) 

ln(l - z 



(3.9) 



ln(l - z 

1 - z 

8tt 2 " 



+ 15 



1 - z 



- 12 



L(i-*) 



6 + 55z + 2(22 - 17 z) ln(l - z) - 4 cJJ In ^~ 



(3.10) 



+ 2 



l-z 



<5(l-z)+45'(l-z)-4 
ln(l - ^ 



16 



1-z 



(1-z) 2 . 
+ 2-6z + 4(l + 3z) ln(l-z). 



(3.11) 



The contributions Ac^, and Ac^'^ are a consequence of the application of the equation 
of motion of the effective theory in the evaluation of the matrix elements of the operators 
given in ( I2.12p . Their explicit expressions read 



Ace- F) = 2c (0) 



'A 2 



A 2 ' 



Ac' 1 -' 



2(l + l„A 4»). 
m 6 7 2 



(3.12) 



We note that the coefficient agrees with the well-known result of [2T]. The coefficients 

and have 
expected relation 



and have been calculated here for the first time. For the first one we confirm the 



dz c, 



(n) 



dz c$™ } , 



n = 0,l,2, ... . 



(3.13) 



Not surprisingly, diverge more strongly than c^ 1 ' at the endpoint. The plus-distributions 
introduced above follow the prescription 



/ dz 


"ln n (l-z)" 


_ (1 - z) m _ 


la 



p (l-^) p 



d p f(z) 
dzP 



(3.14) 



z=l 



where f(z) is an arbitrary test function which is regular at z = 1, and n > 0, m > 1. In 
case the integration does not include the endpoint, we have (c < 1) 



dz 



ln n (l - z) 
(1 - z) m 







(1-z)' 



(3.15) 



We remark that the /i-dependence of Cg has its origin in the MS renormalization of the 



electromagnetic dipole operator. The same is true for the /i-dependence of c£' and c 
It reflects the fact that Ai is not renormalized to all orders in perturbation theory. On the 



,(1,F) 



(1,A) 
A 2 



It 



other hand, for the coefficient of A2 there is an additional //-dependence present in c 
originates in the MS renormalization of A 2 . Indeed, using the renormalization group equation 

7M°», (3-16) 



aisin) (o) (o). 



4tT 



where — — y - \ l—~ - ivvvmvj. mv ^-v lv jn uv lv m v in — w y t 

In the remainder of this section we discuss the numerical impact of the new contributions 
on the total decay rate of B — > X s 7 as well as on its first and second moments. For the O7 
contribution to the total decay rate we find (using numerical values for Cp and C a) 



2Ca [22J , we recover the //-dependence of Acf' A ' 



77\Ej>E 



1 , dT 

dz 

zO 



77 



dz 



77 



1 + 



Ai-9A 2 (» 



2m 2 



Osjp) ( 16 
4tt V 9 _ 



7T~ 



31nA 

m b 



--ln 2 ( 

3 1 

+ ^ 
2m 2 



z ) 



2z - zl) ln(l 



16 



77' 



31nA 

m b 



-In 2 

3 



z ) 



z ) 



-z (30 + 3z - 2zl) 



Az (30 



63^o + 3l2 2 + 5zg) 



4 (30 - 72z + 51^ 2 - 2zl - 3z%) 



+ 27 In 



m b 
2z (54 - 



9(1 -z y 

-f ln 2 (l-^o) + 



49^o + 15^ 2 ) 



ln(l -Zo)> + 



9(1 

A 2 (/i) f 2 



2m 2 



9 



z ) 



17 + 32tt a 



2 (162 - 244z + 113z 2 - 7z%) 
3(1 -zo) 



ln(l - zq) 



Z() 



(3.17) 



The ellipses denote higher order terms in a s and A./m b , and we used z = 2E /m b . I n order 
to get a rough estimate of the size of the power-corrections at 0(a s ) we set /i = m b and use 
the numerical values a s (m b ) = 0.22, m b = 4.6 GeV, A x = -0.4 GeV 2 and A 2 = 0.12 GeV 2 to 
obtain T^lE^i.sGeV /^ri = 0.763 — 0.007 = 0.756 , a —0.9% effect. In the intermediate step 
we have singled out the new 0(a s A 2 /m 2 ) contributions. The effect of the new corrections on 
the rate varies with the cut, from —0.4% at E = to —0.9% at E = 1.8 GeV. For values 
of E > 1.8 GeV the corrections to the rate are significant, a —3% effect for E = 2 GeV, 
however, such high values of E are well outside the range of applicability of the local OPE. 
Moreover, at high E part of the new effect is implicitly contained in the approach of [8]. 

The truncated n-th moment is defined through 



E-y>E 



7/-C/7 



~2~ 



dz z 1 



dT 



77 



dz 



dz 



dT 



77 



' Z / / \J z 

After expanding in a s , the first moment reads 



dz 



(3.18) 



E~,>E 



m b 
2 



Ai + 3A 2 Qu) ajjj) 
2ml 4?r 



46 8, . 

_ + _( 8 _ 9zb + z 3) ln(l-zo) 
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+ ^z (96 - 60z - 22z 2 + 9z 3 ) + 



Ai f 22 (96 - 12(k - 52^ + 25z 3 - 9z*) 



+ 



46 16 (4 - 7z + z, 



27 



^ ln(l - zq) } + t 

9(1 -zq) 1 0) j 2ml 



27(1 - zq) 

A 2 (//)/ 1130 u 4 

- ivv - - 18 In — + - (80 - 33^0 
27 m h 9 



81z 2 + 25^) ln(l - zq) + — z (480 + 42z - 425,? 2 + 8U 3 ) 



+ . . . 



(3.19) 



while the second central moment is given by 



Ai t a s (/i) frn% 
4 



;i-^ ) 2 (l7-2zo-3z 2 )ln(l-z ) 



+ 270 (1 ~ ^ ^ 61 " 898Z ° ~~ 207z2 ° + 144Z °) 



4 



Ai 



27 



(2 - 27^o + 7z 3 ) ln(l - z ) 



+ — (169 + 602 - 780zl - 385zg + 225zj - 54zjj) 



+ 



A 2 (/^) 



9 



(43 - 84^ 



+ 258^ 2 - 292zjj + 75^) H 1 - *b) 



540 



(707 - 25802 + 37502 2 - 13820zg 



+ 145354 - 25924) 



+ ... . 



(3.20) 



Using the same numerical input as above we obtain, for Eq = 1-8 GeV, (Ej) = (m{,/2)(0.978— 
0.001) = 2.246 GeV and (£ 2 ) - (£ 7 ) 2 = 0.0422-0.0036 = 0.0386 GeV 2 , where we have again 
singled out the contribution of 0(a s A 2 /m 2 ) in the intermediate steps. Their effect on the 
truncated first and second central moment is of around —0.1% and —8.5%, respectively. For 



E G [0, 1.8] GeV, the effect on the first moment varies between —0.2% and 
that on the second central moment varies between —3.5% and —8.1 



-0.1%, whereas 



/o. 



It is interesting to compare the 0(a s ) coefficients of Ai in f )3.17p . (I3.19P and ( I3.20p with the 
results of [15] . where the corrections of 0(a s Ai/A 2 ) with A = mb(l — zo) have been computed 
for the cut rate and the first two moments. Expanding our results in A/m^ = 1 — zq and 
keeping only the leading term we reproduce the results of [13]. Fig. [3] summarizes the 
numerical relevance of the terms suppressed by powers of A/rrib in the 0(a s ) coefficients 
of Ai/(2m 2 ) for the same numerical input used above. In the cut rate and in the second 
moment, the leading approximation deviates by roughly +50% and —35%, respectively, 
already at Eq = 2 GeV. In the first moment, the leading approximation is within roughly 
10% of the complete result down to E = 1.8 GeV. In conclusion, the range of applicability of 
the leading order approximation in the expansion in A /mi, is clearly restricted to the region 
E > 2 GeV, where Sudakov logarithms become dominant, see for instance [23] . 

Fig. H] shows the ratios of NLO to leading order coefficients of Ai j2 in the rate and in the 
first two moments as a function of the cut Eq using the same input as above. The NLO 
corrections to A 2 are close to 20%. Note that in the right panel we have not shown a curve 
for the second central moment since A 2 has a vanishing leading order coefficient. 

The second moment represents a powerful constraint on the kinetic expectation value 
/i 2 = — Ai + 0(l/mb). Since the 0(a s ) correction decreases its coefficient by 5 to 9% in the 
range of cuts between and 1.8 GeV, while the 0(a s A 2 /m 2 ) corrections are much smaller, 
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Figure 3: NLO coefficients of Ai (solid curves) and its leading approximation in A/m^ 
(dashed curves) for the decay rate (left panel), the first moment (right panel, lower red 
curves) and second moment (right panel, upper blue curves) as a function of Eq. 




0.0 0.5 1.0 1.5 ' 0.0 0.5 1.0 1.5 



Eq Eq 

Figure 4: Ratio of NLO to leading order coefficients of Ai (left) and A 2 (right) in the rate 

(red solid curves), the first moment (blue dashed curves) and the second moment (black 
dash-dotted curve) as a function of E . 



we expect to extract a higher value of fi^ from radiative moments once the new corrections 
are included. Indeed, using a s at a more appropriate scale of order 1-2 GeV, the extracted \x\ 
gets shifted by approximately +10%. Analogously, the small correction to the first moment 
leads to a roughly 10 MeV positive shift in mj. 

All the above expressions refer to the on-shell scheme for m& and Ai, which is inherent in the 
HQET calculation. In practical applications however one adopts short- distance definitions 
of these parameters, as in the kinetic scheme [23]. In this scheme the new corrections 
are identical, but they additionally induce small 0(a;?/i£ in /m^) contributions. Therefore, 
the above rough estimates hold in the kinetic scheme as well. A complete phenomenological 
analysis of the moments, including all the available contributions, will be presented elsewhere. 
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4 Conclusions 



We have computed the NLO contributions to the Wilson coefficients of dimension five opera- 
tors relevant for inclusive radiative B decays. Our results allow for a more precise evaluation 
of the moments of the photon distribution and will improve the determination of m& and of 
the kinetic expectation value, fi^, from radiative moments. We have estimated that the new 
contributions shift the value of extracted from the radiative moments by approximately 
+10% and that of mb by roughly +10 MeV. The effect on the B — > X s 7 rate is below 1% for 
E < 1.8 GeV. 

We have performed the calculation analytically, using an off-shell matching procedure, 
a method that can be applied to inclusive semileptonic decays as well. The 0{a s ^jm^) 
corrections to the moments of B — > X c iv have been computed numerically [25], however, 
the 0(a s fiQ/ml) corrections are not yet known. We also believe that an analytical result 
might be easier to implement in the fitting codes. 
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